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CLUTTER  MODELS  FOR  SUBSURFACE 
ELECTROMAGNETIC  APPLICATIONS 


David  A.  Hill 

Electromagnetic  Fields  Division 
National  Institute  of  Standards  and  Technology 
Boulder,  CO  80303 


Clutter  models  for  subsurface  electromagnetic  applications  are 
discussed  with  emphasis  on  tunnel  detection  applications. 
Random  medium  models  are  more  versatile  and  require  less 
detailed  information  than  deterministic  models.  The  Born 
approximation  is  used  to  derive  expressions  for  the  incoherent 
field,  and  electric  and  magnetic  dipoles  are  treated  in  detail. 
When  random  inhomogeneities  are  located  in  the  near  field  of  the 
dipole  source,  an  electric  dipole  radiates  a larger  incoherent 
field  than  a magnetic  dipole  because  of  its  larger  reactive 
electric  field. 

Key  words:  clutter;  coherent  field;  electric  dipole;  incoherent 
field;  magnetic  dipole;  random  medium. 


1 . INTRODUCTION 


The  use  of  electromagnetic  (EM)  waves  in  subsurface  applications 
continues  to  increase.  Some  current  applications  are  geophysical 
prospecting  [1-3],  mine  communication  [4,5],  tunnel  detection  [6-8],  and 
mine  detection  [9,10].  The  simplest  earth  model  is  a homogeneous  medium;  a 
horizontally  stratified  medium  [11,12]  is  a better  earth  model  for  some 
applications  [13-16]. 

EM  transmission  measurements  for  tunnel  detection  [7]  and  scattering 
measurements  for  remote  sensing  [17]  indicate  that  the  real  earth  is  often 
inhomogeneous  in  a fairly  complicated  manner.  When  an  EM  wave  propagates 
through  an  inhomogeneous  medium,  it  suffers  attenuation  and  dispersion,  and 
an  incoherent  (clutter)  signal  is  produced  by  scattering.  Such  effects  are 
detrimental  to  tunnel  detection  [7],  but  they  have  not  been  quantitatively 
evaluated.  In  some  applications  of  geophysical  tomography  [18-19],  random 
noise  has  been  added  to  the  received  signal  to  test  the  sensitivity  of  the 
imaging  algorithm.  Such  tests  are  useful,  but  it  would  be  more  realistic  if 
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we  had  mathematical  models  for  clutter  signals  that  represent  the  physics  of 
an  inhomogeneous  earth. 

A typical  borehole- to-borehole  example  for  tunnel  detection  is  shown  in 
figure  1.  The  received  signal  is  made  up  of  a direct-path  signal,  a signal 
scattered  from  the  target  (in  this  case  a tunnel),  and  a clutter  signal 
scattered  from  the  inhomogeneities  in  the  rock.  In  typical  analyses  of  EM 
detection  of  tunnels  [20,21],  the  clutter  signal  is  not  included.  However, 
the  clutter  signal  (sometimes  called  geologic  noise  [22])  is  often  the  main 
limitation  in  tunnel  detection.  In  underground  communication  or  source 
location  [23],  the  target  signal  does  not  exist,  but  the  direct-path  signal 
is  still  contaminated  with  the  clutter  signal.  This  report  is  primarily 
concerned  with  clutter  models  for  use  in  the  analysis  of  tunnel  detection 
and  other  subsurface  applications. 

The  organization  of  this  report  is  as  follows.  Section  2 discusses 
types  of  inhomogeneous  earth  models,  both  deterministic  and  random.  Random 
media  models  offer  more  flexibility,  and  the  mathematical  techniques  for 
analyzing  random  media  are  discussed  in  section  3.  The  analysis  of  an 
electric  dipole  radiating  in  a random  medium  is  contained  in  section  4. 
Numerical  results  are  given  for  both  the  far- field  pattern  and  the  total 
radiated  incoherent  power.  Analogous  results  are  given  for  a magnetic 
dipole  source  in  section  5.  Section  6 contains  an  analysis  for  transmission 
between  electric  dipoles  for  both  cw  and  pulse  excitation.  Conclusions  and 
recommendations  for  further  work  are  contained  in  section  7. 


2.  INHOMOGENEOUS  EARTH  MODELS 

In  this  section,  we  discuss  three  types  of  inhomogeneous  earth  models. 
Deterministic  models  are  useful  when  the  details  of  the  local  earth  medium 
are  well  known,  but  they  require  more  information  than  is  generally 
available  for  most  applications.  Random  media  models  are  more  flexible  and 
require  less  information.  However,  the  EM  fields  propagating  in  random 

media  are  random  variables,  and  they  can  be  described  only  statistically. 

2.1  Deterministic  Models 
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The  most  useful  deterministic  model  is  composed  of  a homogeneous  earth 
with  a single  scatterer  whose  response  can  be  calculated.  In  this  case  the 
scattered  field  is  considered  to  be  the  clutter  signal.  For  example  in  the 
problem  of  surface  location  of  a buried  magnetic  dipole  source,  location 
errors  caused  by  scattering  from  a sphere  [23]  or  a prolate  spheroid  [24] 
have  been  evaluated.  The  geologic  noise  caused  by  a more  complicated 

overburden  scatterer  in  geophysical  prospecting  has  been  evaluated 

numerically  [22]. 

In  these  cases  the  parameters  that  describe  the  scatterer  (dimensions, 
location,  constitutive  parameters,  etc.)  are  deterministic  quantities,  and 
the  scattered  fields  are  deterministic.  Thus  the  scatterer  parameters  are 
chosen  from  a knowledge  of  the  specific  application,  and  each  case  yields 
only  one  deterministic  result  for  the  clutter  signal.  Repeated  calculations 
for  different  parameters  are  required  for  general  clutter  characteristics, 
and  such  repeated  calculations  can  be  tedious.  For  tunnel  detection  where 
the  structure  of  the  rock  is  not  well  known,  we  favor  random  media  models 
because  statistical  properties  of  the  clutter  are  obtained  directly. 
Consequently,  the  remainder  of  this  report  will  deal  with  random  media 
models . 

2.2  Random  Discrete  Scatterers 

Random  media  can  generally  be  grouped  into  three  categories  [25]: 
random  scatterers,  random  continua,  and  rough  surfaces.  We  will  not 
consider  rough  surfaces  in  this  report  because  we  assume  that  volume  clutter 
is  more  important  than  surface  clutter  in  most  tunnel  detection 

applications.  However,  for  cases  where  the  source  or  receiver  is  located 
near  the  earth  surface  or  near  a subsurface  layer  with  a rough  boundary, 
rough  surface  clutter  could  be  important. 

Typical  examples  of  discrete  scatterers  are  rain  [26]  and  hail. 
Discrete  scatterers  can  be  characterized  by  their  bistatic  radar  cross 

A A A 

section  cr^^(o,i),  where  i is  a unit  vector  in  the  direction  of  propagation 

A 

of  the  incident  wave  and  o is  the  direction  of  the  scattered  wave  [25,  Sec. 
2-1].  (Vectors  and  dyads  are  indicated  by  boldface  throughout  this  report.) 
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For  a random,  uniform  distribution  of  identical  scatterers,  the  other 
important  parameter  is  the  number  p of  scatterers  per  unit  volume. 

For  a tenuous  distribution  of  scatterers  in  the  far  field  of  the 
source,  such  as  a weather  radar  [26],  the  single  - scattering  approximation 
yields  a fairly  simple  volume  integral  expression  for  the  scattered  power 
[25,  Sec.  4-1].  The  cases  of  a dense  distribution  or  a near-field  source 

are  more  complicated. 

2.3  A Random  Continuum 


In  this  report  we  define  a random  continuum  as  a medium  whose 
permittivity  is  a random  continuous  function  of  position.  In  some  cases, 
such  as  atmospheric  turbulence  [27],  the  permittivity  is  also  considered  to 
be  a function  of  time.  In  modelling  the  earth,  we  assume  that  the 
permittivity  e(r)  is  a random  function  of  position  r only. 

We  write  the  earth  permittivity  as  a sum  of  the  average  value  and  a 

fluctuation : 

e (r)  = e Q [ 1 + e^r)  ] , (1) 


where  is  the  fluctuation  with  zero  average 

<e1 (r)>  = 0 (2) 

and  < > represents  the  ensemble  average  [28].  For  lossless  media,  and 

are  real,  but  for  lossy  media  we  make  them  complex.  If  the  medium  is 

statistically  homogeneous  and  isotropic,  then  the  covariance  <e^(r^) 

is  a function  only  of  the  difference  r , = |r,|  = I r.  - r~ I : 

d d 1 z 

i (ri ) e^(r2)>  = (rd) > (3) 

where  superscript  * denotes  complex  conjugate. 
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The  equivalence  between  continuous  and  discrete  scattering  models  has 

been  studied  by  Besieris  [29].  As  we  will  see  in  section  3,  the  scattering 

cross  section  per  unit  volume  can  be  written  in  terms  of  the  Fourier 

transform  of  B . Various  analytical  forms  of  B (r,)  have  been  used  [25,  Ch. 
e e d 

16],  and  anisotropic  forms  of  B^(r^)  have  been  used  for  anisotropic  media 

such  as  vegetation  [30,31].  In  this  report  we  use  only  the  isotropic  form 
B^(r^)  because  it  is  simpler,  but  anisotropic  forms  might  be  more  realistic 

for  modelling  nonspherical  inhomogeneities  or  fractures. 


3.  MATHEMATICAL  METHODS  FOR  RANDOM  MEDIA 

An  EM  wave  propagating  in  a random  medium  can  be  expressed  as  the  sum 
of  an  average  field  and  a fluctuating  field.  The  average  field  is  also 
called  the  coherent  field,  and  the  fluctuating  field  is  also  called  the 
incoherent  field  [25,  Sec.  4-4].  For  example,  we  can  write  a scalar 

component  of  the  total  time -harmonic  electric  field  E^(r)  as 

Eu(r)  = <Eu(r)>  + Efu(r) , (4) 

where  <E^(r)>  is  the  average  (coherent)  field  and  E^  (r)  is  the  fluctuating 

(incoherent)  field.  The  time  dependence  exp(-iwt)  is  suppressed  throughout 
this  report. 

3.1  Coherent  Field 

When  an  EM  wave  propagates  through  a random  medium,  the  average 

(coherent)  field  is  generally  attenuated  and  slowed  by  scattering.  However, 
if  the  randomness  is  very  weak,  then  the  coherent  field  can  be  approximated 
by  the  incident  field  that  would  propagate  in  a uniform  medium.  The 

approximation  of  the  coherent  field  by  the  incident  field  is  known  as  the 

Born  approximation  [25,  Sec.  4-1]. 
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Karal  and  Keller  [32,33]  have  used  a perturbation  technique  to  study 
propagation  in  a random  continuum.  They  have  derived  a dispersion  equation 
for  the  effective  wavenumber  for  plane -wave  propagation  and  have  derived 
expressions  for  the  effective  dielectric  constant  and  conductivity.  Roth 
and  Elachi  [34]  have  used  the  dispersion  equation  of  Karal  and  Keller  to 
generate  numerical  results  for  loss  tangent  as  a function  of  the  parameters 
of  a random  continuum  with  an  exponential  correlation  function. 

When  the  medium  has  strong  fluctuations,  perturbation  methods  are  no 
longer  valid.  Dyson's  equation  for  the  coherent  field  [34]  is  exact  in 
principle,  but  various  approximate  methods  [34-39]  are  required  to  obtain  a 
solution . 


3.2  Incoherent  Field 

By  definition  the  average  value  of  the  incoherent  field  is  zero.  For 
example,  from  (4)  we  can  see  that  <E^>  = 0.  However,  if  we  take  the 

average  of  the  square  of  the  absolute  value  of  (4) , then  we  can  write 

< | E |2>  - | <E  >|2  + <|E„  |2>.  (5) 

1 u 1 u 1 f u 1 

2 

If  we  use  Ishimaru's  terminology  [25,  Sec.  4-4],  then  < | | > is  the  average 

2 2 
intensity  <I>,  | <E  >|  is  the  coherent  intensity  1^,  and  <|E^J>  is  the 

incoherent  intensity  I..  Thus  (5)  could  be  written 


<I>  = I + I. . (6) 

c 1 

Ishimaru's  definitions  are  actually  applied  to  scalar  fields,  but  we  can  use 
the  same  terms  for  scalar  components  E^  of  the  electric  field  or  the  vector 

2 

field  E.  If  we  want  the  intensities  to  have  units  of  W/m  , then  each  term 
in  (5)  or  (6)  should  be  multiplied  by  the  intrinsic  impedance  rj  of  the 
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medium. 


In  most  cases  we  will  be  interested  in  the  ratio  I. /I  , and  the 

l c 

intrinsic  impedance  will  cancel  out. 

In  a weakly  random  medium,  the  coherent  intensity  is  nearly  the  same  as 
the  intensity  in  absence  of  randomness,  and  the  incoherent  intensity  1^  can 

be  considered  the  intensity  of  the  clutter  signal.  Thus  the  main  problem  in 
clutter  applications  is  the  evaluation  of  the  incoherent  intensity  I = 

2 

<|Eful  >•  The  simplest  methods  for  calculating  I ^ use  the  Born 

approximation  [25,  Sec.  2-6],  and  we  will  use  the  Born  approximation 
throughout  the  remainder  of  this  report.  There  are  more  sophisticated 

methods,  such  as  the  Bethe - Salpeter  equation  [25],  for  calculating  1^,  but 

the  earth  medium  has  not  been  well  characterized  in  tunnel  detection,  and 
there  is  no  point  in  making  the  initial  calculations  overly  complicated. 

We  can  start  the  analysis  by  writing  the  source -free  wave  equation  for 
the  electric  field  E, 

V X V X E - 1 + €1)E  = 0,  (7) 

where  \i ^ is  the  permeability  of  the  medium  and  the  permittivity  has  been 
written  in  the  form  shown  in  (1) . We  now  write  (7)  in  the  equivalent  form, 

V X V X E - k^E  = k^e  E,  (8) 

2 2 

where  k^  = u>  Mq€q-  We  can  write  the  total  field  E as  the  sum  of  the 
incident  field  E^  and  the  scattered  field  E^: 

E = EQ  + Er  (9) 

The  incident  field  satisfies 
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V x V x Eq 


k2E 
0 0 


0. 


(10) 


If  we  subtract  (10)  from  (8),  we  obtain 


V x V x E^ 


koEi  - k0£iE- 


(ID 


We  assume  that  we  know  the  incident  field  which  satisfies  (10) . The 
scattered  field  E^  satisfies  the  wave  equation  with  a source  term  on  the 

right  side  as  given  by  (11).  The  problem  is  that  the  source  term  in  (11)  is 
not  known  because  E is  not  known.  In  the  Born  approximation  [25,  2-6],  we 
approximate  E by  the  incident  field  E^  in  the  source  term  in  (11).  Thus  the 

scattered  field  approximately  satisfies  the  inhomogeneous  wave  equation, 


V X V X E^ 


koEi  koeiEo 


(12) 


Since  the  right  side  of  (12)  is  known,  we  can  write  E^  in  integral  form 
using  the  Green's  function  approach  [40] 


Vr)  - k0 


J Mr') 

V 


G (r , r ' ) 


Eo(r'} 


dr' 


(13) 


where  G(r,r')  is  the  dyadic  Green's  function  and  the  integration  volume  V is 
all  space.  In  some  remote  sensing  cases,  G(r,r')  could  be  the  dyadic 
Green's  function  for  a layered  medium  [29,30],  but  we  consider  only  a 
homogeneous  medium  (constant  k^)  for  the  deterministic  portion  of  the 

problem.  In  this  case,  G(r,r')  is  [40] 


G(r , r ' ) 


zi  (l  + h vv)  *• 

ko 


(14) 
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where  if>  = expCik^lr  - r'|)/|r  - r'|  and  1 is  the  unit  dyad.  For  points 

where  r = r',  either  G or  the  integral  in  (13)  must  be  appropriately 
modified  [41,42].  We  can  also  write  (13)  in  component  form  for  a scalar 
component  E^  of  E^: 


Elu(r)  - k0 


I 

V 


ex(r') 


G (r , r ' ) 
uv 


E0v(r>) 


dr' 


(15) 


where  G is  a scalar  component  of  G and  the  repeated  subscript  v indicates 

summation  over  the  three  values  of  v.  From  (15)  we  can  see  that  <E-  > = 0 

lu 

because  <£^>  = 0. 

2 

To  evaluate  the  incoherent  intensity,  we  need  to  evaluate  <|E^  | If 

we  multiply  (15)  by  its  complex  conjugate  and  take  the  average  value,  we 
obtain 


<|£  (r)|2>  - |k  I4  J ; B ( | r ' 

v V 


r"l)  Guv(r,r')E0v(r’). 


G (r  , r ' 
uw 


)E0w(r' 


) dr'  dr' 


(16) 


where  we  have  used  (3)  to  introduce  into  (16) . In  the  modified  Born 

approximation  [43,44],  the  Green's  function  in  (16)  is  replaced  by  the 
effective  Green's  function  for  the  random  medium. 

Generally  the  double  volume  integral  in  (16)  is  difficult  to  evaluate, 
but  when  the  random  volume  is  in  the  far  field  of  the  source  and  observer, 
the  integrals  simplify  considerably  [25,27].  For  the  geometry  in  figure  2, 
Ishimaru  [25,  Sec.  16.1]  writes  the  received  power  P^as 

A2G  (i)G  (o)  A a 

P “ U A A A <7(0,  i)  dV,  (17) 

V (4tOZR^R2 
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where  P^_  is  the  transmitted  power,  and  are  the  gains  of  the 

transmitting  and  receiving  antennas,  and  are  the  distances  shown  in 

figure  2,  V is  the  volume  shown  in  figure  2,  and  A is  the  wavelength.  We 
assume  here  that  the  medium  is  lossless,  but  an  exponential  loss  factor  can 
be  included  when  necessary.  The  bistatic  scattering  cross  section  per  unit 

A A 

volume  a(o,i)  can  be  derived  in  terms  of  B^  in  integral  form  [25,  16.2] 

<7(0, i)  - *kQ  sin2*  4>£(ks)/2,  (18) 


where  $ (K)  = r f B (r,)  exp(iK»r,)  dr,, 

e . „ . 3 J e d d d 

( An)  °o 

= 2k^  sin(0/2), 

A A 

and  the  unit  vectors  o and  i are  shown  in  figure  2.  The  scattering  angle  9 
and  the  polarization  angle  \ are  shown  in  figure  3. 

There  are  many  possible  forms  for  B , but  the  exponential  form  is  one 

which  is  mathematically  convenient  [45] 

B (r  ) = <e?>  exp ( - r /a) . (19) 

e d 1 d 

2 

This  model  is  characterized  by  the  variance  <e^>  and  the  correlation 
distance  a.  The  integral  in  (18)  is  easily  evaluated  to  yield 

4 3 2 2 

kna  sin  y <e  > 

y1 T . (20) 

2tt(1  + 4k^a  sin (0/2)) 

When  the  correlation  distance  a is  small  compared  to  wavelength  A,  the 
cross  section  a reduces  to 
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(21) 


o{6)  ~ k^a3  sin2*  <e^>/(2n) . 

The  expression  in  (21)  represents  Rayleigh  scattering  where  the  scattered 
power  is  proportional  to  frequency  to  the  fourth  power  times  the  correlation 
length  to  the  third  power. 


4.  ELECTRIC  DIPOLE  RADIATION 

In  the  previous  section  we  noted  that  the  double  volume  integral  in 
(16)  is  difficult  to  evaluate  except  when  the  random  medium  is  in  the  far 
field  of  both  the  source  and  receiver.  Dorfman  [46]  has  treated  the  case 
where  an  electric  dipole  is  located  within  a spherical  random  medium,  and  he 
was  able  to  evaluate  the  intensity  in  the  far  field.  In  this  section  we 
elaborate  on  his  treatment,  and  in  section  5 we  treat  the  analogous  case  for 
a magnetic  dipole  source. 

4.1  Far-field  Pattern 

The  geometry  for  a z-directed,  electric- dipole  source  of  moment  IL  is 
shown  in  figure  4.  The  randomness  of  the  medium  is  confined  to  a spherical 
shell,  b < r < c,  where  the  permittivity  is  given  by  (1).  The  spherical 
shell  model  is  idealized,  but  it  allows  us  to  examine  near- field  effects  for 
the  transmitting  dipole.  Following  Dorfman  we  assume  that  the  medium  is 
lossless,  but  loss  could  be  added  to  the  model.  Since  we  use  the  Born 
approximation,  the  coherent  far  field  has  only  a 6 component  of  electric 
field  En  given  by  [47] 

U o 

- iw/^  IL  sin# 

E06  ~ W exp(ikQr)  . (22) 

The  coherent  intensity  I is 

y c 
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(23) 


(w/i qIL  sin# ) 

k - 'Eod  - — — T 

(47rr ) 


Following  Dorfman  [46]  we  approximate  B in  the  form, 


Bf(rd)  - <£2>  P «(rd) , 


(24) 


where  8 is  a Dirac  delta  function  and  2 plays  the  role  of  a correlation 
length.  The  form  of  is  physically  realistic  only  when  the  incident  field 

is  nearly  constant  over  the  distance  i.  The  main  advantage  of  (24)  is  that 

it  allows  the  integral  in  (16)  to  be  evaluated  analytically.  If  the  actual 

correlation  function  is  r(r,),  then  i is  obtained  from 

d 


1 - S r(rd)  drd- 


(25) 


3 

For  example,  if  T is  an  exponential  function  exp(-r^/a),  then  f is 


3 r r 2 3 

i = J exp(-r^/a)  dr^  = 47r  J exp(-r^/a)  r^  dr^  = 87ra  . 


0 


(26) 


If  we  substitute  (24)  into  (16)  and  carry  out  the  r''  integration,  we 
obtain 


< 


Elu(r)|2> 


= k 


4 


«2i> 


I3  f | G (r,r')E.  (r') |2 
^ 1 uv  Ov 


dr ' 


(27) 


where  V is  the  spherical  shell  volume.  When  the  observation  point  is  in  the 
far  field  (large  k^r)  , the  Green's  function  in  (14)  or  (27)  simplifies 

considerably,  and  (27)  simplifies  to 
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(28) 


< I E]_0  (r)  I 2> 

<t> 


2 3 2 

VZ&<e^> 

(47rr ) ^ 


f \* 
V * 


EQ(r')|  dr', 


A A ^ 

where  5 and  <f>  are  unit  vectors.  In  the  far  field,  <|'E^  (r)  | > ~ 0. 
For  an  electric  dipole  source,  E^(r')  is  [47] 


Eo(r<)  - 


-iw/z^IL  exp(ikQr')  A _2 


47rr ' 


lr'  [T 


+ 


ik  r'  2 

0 (ikQr') 


] cos 6 ' 


(29) 


+ O'  [1  - tT—  + — 1 T]  sin0'}. 

lk0r  fik0r’) 


If  we  substitute  (29)  into  (28) , we  can  carry  out  the  volume  integration 
analytically.  The  details  are  included  in  Appendix  A,  and  the  result  is 

<|E1S |2>  - k3i3  <e]>  E2  , (30) 

<t>  <f> 

where  F = [k  (c  - b)(cos^/9  + 8 sin^$)  + k/(b  ^ - c ^)(3  cos^d  + 4 sin^0) 
v U U 

+ k"  3 (b"  3 - c"3)  (23  cos  2 $ + 4 sin20)]/(6O  ir)  , 

F - [kQ(c  - b)  + Sk'bb"1  - c'1)  + 15kQ3(b'3  - c'3)]/(60  k)  , 
and  Eq  = IL)/(47rr).  The  incoherent  intensity  I ^ is 

h - <iei^i2>  + <i  Ei/>  - klp  <£i>  Eo  (h  + v-  (31) 
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The  ratio  of  the  incoherent  to  coherent  intensities  is  also  of  interest 
because  it  gives  a measure  of  the  degradation  of  the  transmitted  signal 
between  a pair  of  electric  dipoles: 

I. /I  = 3 <e3>  (F.  + F ) /sin2# . (32) 

1 c 0 1 0 <f> 

The  results  in  (31)  and  (32)  represent  Rayleigh  scattering  where  the 

incoherent  intensity  is  proportional  to  third  power  of  the  correlation 

3 4 

length,  1 . This  is  similar  to  (21).  However,  the  dependence  is 

obtained  only  when  the  k c term  is  dominant  in  F and  F , . This  is  as 

U u (p 

expected  because  this  is  the  case  where  far-field  scattering  is  dominant. 
When  near- field  scattering  is  important,  then  the  k^  dependence  is  to  some 

-3  -3 

lower  power.  When  the  quasi-static  term,  kg  b , is  dominant,  the 

scattering  is  independent  of  k^ . If  we  let  b approach  zero,  then  becomes 

infinite.  This  physically  unrealistic  result  is  a result  of  using  the  delta 
correlation  function  in  (24).  If  we  had  chosen  a finite  correlation 
function,  then  I would  have  remained  finite  as  b approached  zero  [46,48]. 

However,  it  is  difficult  to  evaluate  the  volume  integral  for  finite 
correlation  functions.  The  results  in  this  section  are  still  valid  as  long 
as  b is  somewhat  greater  than  i. 

Numerical  results  for  the  incoherent  intensity  are  shown  in  figures  5 

2 

and  6.  The  results  are  normalized  to  CE^  to  make  them  dimensionless,  and 

3 3 2 

the  normalization  constant  C = k^i  <e^>.  In  figure  5,  far-field  scattering 

is  dominant,  and  the  maximum  occurs  in  the  broadside  direction  (8  — 90°)  as 
it  does  for  the  coherent  component.  In  figure  6,  near-field  scattering  is 
dominant  because  of  the  small  value  of  k^b . Here  the  maximum  occurs  off  the 

end  of  the  dipole  ( 6 = 0°).  An  interesting  feature  of  both  cases  is  that 
the  <f>  component  is  independent  of  8.  The  cross -polarized  component  is 
normally  zero  for  far-field,  plane-wave  scattering  in  the  Born  approximation 
[49],  but  not  for  the  near-field  case  considered  here. 
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4.2  Total  Radiated  Power 


The  total  radiated  power  can  be  obtained  by  integrating  the  intensity 
over  a sphere  at  infinity.  The  coherent  power  P^  is 


po  - 

27T  7T 

f f I r sin#  d#  d<fi . (33) 

Zrjo  o o c 

If  we  substitute  (23)  into  (33)  and  carry  out  the  integrations,  we  obtain 


po  - 

r)0(k0IL)2/(12*).  (34) 

The  incoherent  power  P^  is 


P1  - 

2n  n 

> J J I.  r sin#  d#  d<f> . (35) 

Zr)0  0 o 1 

If  we  substitute  (30)  and  (31)  into  (35)  and  carry  out  the  integrations,  we 


obtain 

pi  - 

WL>2  & <•?>  [.l  . vi.ri  -i 

9 [5k  (c  - b)  + 5k  (b  - c ) 

360  n 

(36) 

-3  -3  -3 

+ 19kg  (b  - c )]. 

The  results  in  (31)  and  (36)  are  similar  to  those  of  Dorfman  [46],  but  the 
coefficients  are  not  the  same.  It  is  not  possible  to  check  Dorfman' s 
derivation  because  he  does  not  give  any  intermediate  steps. 

The  ratio  of  incoherent  to  coherent  power  can  be  regarded  roughly  as  a 
clutter- to  - signal  ratio  and  is  given  by 
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0 


1 3,3  2 

V <€1> 

30  7T 


1„  -1 


[5kQ(c  - b)  + 5kQ^(b 


-1  -3  -3 

c ) + 19k  (b 


- c 


)].  (37) 


Because  we  have  assumed  that  the  observation  point  is  in  the  far  field  of 
both  the  source  and  the  random  region,  P^/P^  is  independent  of  r.  Normally 

the  ratio  of  the  incoherent  to  coherent  intensity  would  increase  with  r if 
the  entire  medium  were  random  [25,  Ch.  6]. 


5.  MAGNETIC  DIPOLE  RADIATION 

In  this  section  we  extend  the  results  of  the  previous  section  to  a 
magnetic  dipole  source.  A magnetic  dipole  source  (small  loop)  has  been 

suggested  for  the  detection  of  long  conductors  in  tunnels  [50], 

5.1  Far-field  pattern 

The  geometry  for  a magnetic  - dipole  source  is  the  same  as  that  in  figure 

4 except  that  the  electric-dipole  source  is  replaced  by  a z-directed 

magnetic-dipole  source  of  moment  IA.  Here  we  will  compute  the  magnetic 

field  because  we  anticipate  reception  with  a magnetic  dipole  in  borehole-to- 

borehole  transmission.  The  coherent  far  field  has  only  a 6 component  of 

magnetic  field  H_  given  by  [47] 

U v 


IA  sin# 
4?rr 


exp ( ikQr ) . 


The  coherent  intensity  I is 

c 


I 

c 


(k^IA  sin0)^ 
(47rr ) ^ 


(38) 


(39) 
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If  we  again  use  the  Born  approximation  and  the  Green's  function 
approach  [40] , we  can  derive  the  magnetic  field  analog  to  (13)  for  the 
incoherent  field  H^: 


-iwe 

Hpr)  - — — f c^r')  W x EQ(r')  dr'  , 


(40) 


where  V*  = expCik^lr  - r'|)/|r  - r'|.  As  in  the  previous  section,  we 

approximate  B^  by  (24)  and  carry  out  the  delta  function  integration.  In 

addition  we  make  the  far-field  approximation  and  obtain  the  magnetic  - field 
analog  to  (28) : 


3 3 2 

9 coe  k i <e,>  a a 9 

<|H  (r)|>  = f \6  • [r  x E0(r')]  | dr 

4 (W)  V l 


In  the  far  field,  <|H^  (r)|  > ~ 0. 


For  a magnetic  dipole  source,  E^(r')  is  [47] 


(41) 


EQ(r')  = 4> 


w/i0kQIA  exp(ikQr') 
47rr ' 


(1  + sin0  ' . 


(42) 


If  we  substitute  (42)  into  (41) , we  can  carry  out  the  volume  integration 
analytically.  The  details  are  included  in  Appendix  B,  and  the  result  is 

<|H19|2>  = k^i3  H2  Ce,  (43) 

4>  4> 

where  Gg  - [kQ(c  - b)  + - c" 1) } /(12tt)  , 

G,  = cosh  [kA(c  - b)  + kh(b  2 - c 
(p  DU 
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and  = (w6Qk^IA)/(47rr)  . The  incoherent  intensity  I is 

h - <ihi42>  + <iy 2>  - koi3  <£h  Ho  y + y • <44) 

The  ratio  of  the  incoherent  to  coherent  intensities  is  also  of  interest 
because  it  gives  a measure  of  the  degradation  of  the  transmitted  signal 
between  a pair  of  magnetic  dipoles: 

I. /I  = kV  <e*>  (G.  + GJ/sinV  (45) 

i c 0 10  (f>  ' 

The  main  difference  between  the  magnetic  and  electric  dipole  results  is 

-3  -3 

that  the  magnetic  dipole  result  in  (43)  contains  no  quasi-static  b 

term.  This  is  so  because  the  primary  electric  field  of  the  magnetic  dipole 
/ -3 

in  (42)  contains  no  quasi-static  r'  term.  Thus  the  near-field  scattering 
from  random  irregularities  near  the  source  is  much  weaker  for  a magnetic 
dipole  source,  and  this  is  a potential  advantage  for  magnetic  dipole  (small 
loop)  antennas. 

Numerical  results  for  the  incoherent  intensity  are  shown  in  figure  7. 
The  maximum  scattering  occurs  off  the  end  of  the  magnetic  dipole  (0  = 0° , 
90°).  Again  we  see  that  there  is  a cross -polar ized  <f>  component. 

5.2  Total  Radiated  Power 


The  total  radiated  power  can  be  obtained  by  integrating  the  intensity 
over  a sphere  at  infinity.  The  coherent  power  is 


P = 

0 2 


ij  2n  n 

— | J I r sin0  d9  d <f> 

0 0° 


(46) 


If  we  substitute  (39)  into  (46)  and  carry  out  the  integrations,  we  obtain 
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(47) 


PQ  - (k^IA)  2/(  127r)  . 


The  incoherent  power  P^  is 


r?  2tt  7T 

P t - f-  f f I.  r sin#  d#  &4>.  (48) 

0 0 


If  we  substitute  (43)  and  (44)  into  (48)  and  carry  out  the  integrations,  we 
obtain 


n2TAN2  ,3„3  2 

r?  (k  IA)  k I <ei>  -1  _i  _l 

Px  - -y— ^ — [k0(c  - b)  + k0L(b  - C 1)]. 

72  7T 


(49) 


The  ratio  of  incoherent  to  coherent  power  can  be  regarded  as  a clutter- 
to-signal  ratio  and  is  given  by 


2 3 2 

P k i <ei>  -i  -i  _i 

= ^ - b>  +k0  (b  - ^ )]■ 


(50) 


The  power  ratio  in  (50)  agrees  with  the  first  two  terms  of  the  electric 
dipole  result  in  (37)  . The  difference  is  that  the  electric  dipole  result  in 

-3  -3 

(37)  contains  an  additional  quasi-static  term  proportional  to  k^  (b 
-3 

c ) . Since  the  magnetic  dipole  result  in  (48)  does  not  contain  this  term, 
a magnetic  dipole  does  not  scatter  as  much  power  when  near-field  scattering 
is  dominant  (small  k^b) . Numerical  results  for  the  power  ratios  are  shown 

in  figure  8 for  both  a vertical  magnetic  dipole  (VMD)  and  a vertical 
electric  dipole  (VED) . For  large  k^b , the  results  are  the  same  for  the  two 

dipole  types.  For  small  k^b , the  incoherent  power  rises  rapidly  for  the  VED 

because  of  the  additional  quasi-static  term  in  (37). 

If  we  write  the  two  power  ratios  for  the  small  k^b  limit,  we  obtain 
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(VMD) , 


2 3 2 

pi  V <£i> 


0 


67T 


2 

19  <e^> 
30  7T 


(VED) . 


(51) 


It  is  clear  from  (51)  that  the  incoherent  power  increases  much  faster  (b 
dependence)  for  decreasing  b for  the  VED  than  for  the  VMD.  Also,  the  power 
ratio  is  independent  of  frequency  for  the  VED,  but  is  proportional  to 
frequency  squared  for  the  VMD.  Both  of  these  comparisons  indicate  a 

possible  advantage  for  the  VMD  when  the  source  is  located  within  an 
inhomogeneous  region.  (By  reciprocity  the  same  conclusions  apply  when  the 
receiving  antenna  is  located  within  an  inhomogeneous  region.)  It  is 
interesting  that  similar  conclusions  hold  for  near-field  conduction  losses. 
Input  impedance  calculations  for  a VMD  within  a spherical  [51]  or 

cylindrical  [52]  cavity  in  a conducting  medium  show  the  same  b ^ dependence 
as  (51) . Input  impedance  calculations  for  a VED  close  to  a conducting  half 

-3 

space  [53]  show  the  same  b dependence  as  (49).  The  frequency  dependences 
for  dipoles  in  a lossy  medium  are  somewhat  different  from  those  in  (51) , but 
even  there  the  VED  shows  more  loss  at  low  frequencies  [52]. 


6.  TRANSMISSION  BETWEEN  VERTICAL  DIPOLES 

If  the  transmitting  and  receiving  antennas  are  located  in  an  infinite 
random  medium,  then  the  integrals  arising  from  the  Born  approximation  are 
difficult  to  evaluate  and  might  not  even  converge.  However,  by  neglecting 
backscatter  and  by  making  far-field  approximations,  Ishimaru  [25,  Ch.  6]  has 
developed  some  useful  approximations.  In  this  section  we  specialize  his 
results  to  two  vertical  dipoles  in  a random  continuum.  In  section  6.1  we 
consider  cw  transmission,  and  in  section  6.2  we  consider  pulse  transmission. 

6.1  CW  Transmission 
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The  geometry  for  transmission  between  vertical  dipoles  is  shown  in 
figure  9.  Ishimaru  writes  the  ratio  of  the  received  coherent  power  to 

the  transmitted  power  P^_  as 


P 

c 

P 

t 


.2  GG 

A t r -7 

2 2 6 

(4tt)  L 


(52) 


where  G and  G^_  are  the  gains  of  the  transmitting  and  receiving  antennas, 

and  e represents  attenuation  of  the  medium  due  to  scattering.  For  a lossy 

medium,  e would  also  include  absorption  loss.  For  the  geometry  in  figure 
9,  we  have  L = d/cosa  and  G = G^  = (3/2)cosa.  The  results  in  this  section 

apply  to  both  electric  and  magnetic  dipoles  because  both  have  the  same  far- 
field  pattern,  and  we  can  rewrite  (52)  as 


P 

c 

P 

t 


2 2 4 

A (3/2)  cos  o -7 

2 2 e 
(4tt)  d 


(53) 


An  approximate  ratio  of  the  incoherent  power  P^  to  the  coherent  power 
P^  is  given  by  Ishimaru: 


P./Pc  = 7 - poh, 


(54) 


where  the  total  scattering  cross  section  per  unit  volume  per ^ is  given  by  the 
integral  of  the  bistatic  cross  section  over  47r  steradians  [25] 


A A 

PO  = f °(°’  *■)  dQ> 

47T 


(55) 


and  o is  given  by  (18) . For  the  exponential  model  of  given  in  (19) , o is 
given  by  (20).  In  general  the  integration  in  (55)  is- difficult  to  carry 
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out,  but  approximate  results  for  the  limits  of  small  kna  and  large  kna  are 

u u 

given  in  Appendix  C.  If  these  results  are  substituted  into  (54),  the 
following  results  are  obtained  for  the  ratio  of  the  incoherent  to  coherent 
power : 

p (4/3)<ep>(k0a)3kQL,  kQa  « 1 

^ » ( (56) 

C (l/2)<£3>(k0a)k0L,  kQa  » 1. 

The  approximate  results  in  (56)  also  agree  with  Ishimaru's  plane-wave 
results,  and  for  this  weak  scattering  case  the  power  ratio  is  proportional 
to  the  path  length.  As  usual  the  power  ratio  is  also  proportional  to  the 
2 

variance  <e^>. 

6.2  Pulse  Transmission 

Pulse  transmission  is  also  of  interest  in  tunnel  detection,  and  arrival 
time  and  pulse  dispersion  are  useful  quantities  in  tomographic  imaging  [7]. 
Ishimaru  [25,  Sec.  6.5]  has  obtained  an  approximate  result  for  pulse 
propagation  between  a pair  of  antennas,  and  we  apply  his  result  to  the 
geometry  in  figure  9. 

The  transmitted  pulse  has  energy  and  varies  with  time  as  the  delta 

function  <$(t).  The  received  coherent  intensity  I (t)  has  the  same  time 
variation  with  a time  delay: 

h(t)  - Pc  E0  5(t  - (57) 

where  c^  is  the  velocity  in  the  background  homogeneous  medium  and  P^  is 
given  by  (53) . 

The  received  incoherent  intensity  I^(t)  is  obtained  approximately  by  a 
saddle  point  evaluation  [25,  Sec.  6.5]: 
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(58) 


I.(t)  - 2(,O1/2PcE0Vasc0 


exp [ - ( cq/L) ( 8a^  - pa  ^ L)(t 


L/c 0)] 


The  result  in  (58)  applies  for  t > L/c^,  and  I^(t)  = 0 for  t < L/c^.  Lor 
large  particle  size  D,  a and  the  exponential  time  constant  r are  given  by 


a 

P 


ttD  2 

(2A)  and  T 


L/(8c0a  ) 


2 2 2 
LA  /(2tt  c D ) . 


(59) 


The  total  incoherent  energy  E. 


in  the  received  pulse  is  given  by 


Ei  = f Ii(t)  dt  « PcE0pasL  (jt1/2/4).  (60) 

The  ratio  of  the  incoherent  to  coherent  energy  in  the  received  pulse  is 

E./Ec  = (tt1/2/4)  pas L.  (61) 

This  result  is  the  same  as  the  cw  result  in  (54)  except  for  a constant 
factor  / 4)  . 

For  small  particle  sizes  (D  « A) , and  the  exponential  time  constant 
are  given  by 

a » ' 2,77  « 1.12  and  r « L/(8c  a ) * L/(9c  ) . (62) 

P (n/2)Z  U P U 

In  obtaining  (62)  from  Ishimaru's  result,  we  have  made  the  assumption  that 
the  beamwidth  for  small  particle  scattering  is  tt/2 . This  result  for  the 
decay  time  constant  is  quite  different  from  the  result  in  (59) , but  in  both 
cases  the  length  of  the  incoherent  pulse  is  proportional  to  the  path  length 
L.  Thus  for  either  large  or  small  particles,  the  relative  size  of  the 
incoherent  pulse  grows  in  both  amplitude  and  length  as  the  path  length  L 


23 


increases.  This  is  in  agreement  with  the  cw  description  of  the  incoherent 
wave  in  Ishimaru  [25,  Ch . 6], 

7.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  represents  a first  effort  at  modelling  clutter  in 
subsurface  electromagnetic  applications,  particularly  for  tunnel  detection. 
While  deterministic  earth  models  [22]  can  be  used  for  modelling  clutter,  it 
appears  that  random  media  models  [25]  are  more  versatile  and  require  less 
information  about  the  inhomogeneous  earth.  The  Born  approximation  is  the 
simplest  approach  to  computing  the  incoherent  field  that  is  generated  in  a 
random  medium,  but  even  with  the  Born  approximation  the  calculations  are 
difficult  for  real  antenna  applications  where  the  fields  are  not  plane 
waves.  The  main  difficulty  arises  in  evaluating  the  volume  integrals  of  the 
type  discussed  in  Section  3.2. 

However,  useful  results  can  be  obtained  when  only  the  transmitting 
antenna  is  located  close  to  the  random  medium,  and  specific  results  for 
electric  and  magnetic  dipoles  are  derived  in  Sections  4 and  5.  The  pattern 
of  the  incoherent  intensity  radiated  by  either  type  of  dipole  is  different 
from  that  of  the  coherent  intensity.  A practical  result  is  that  electric 
dipoles  are  more  strongly  affected  by  near-field  inhomogeneities  than  are 
magnetic  dipoles.  Similar  results  have  been  previously  obtained  for  the 
effects  of  near-field  conduction  losses  for  electric  and  magnetic  dipole 
sources  [ 52 ] . 

There  are  many  unsolved  problems  in  clutter  analysis.  The  isotropic 
correlation  functions  that  we  have  used  are  not  applicable  to  flat  or 
fracture  irregularities,  and  fractures  are  thought  to  be  the  most  important 
scatterers  in  tunnel  detection  applications  [54],  It  remains  to  be  seen 
whether  an  anisotropic  correlation  function  or  a totally  new  approach  is 
required  to  model  fractures.  Even  wl  n clutter  is  adequately  modelled,  the 
question  of  how  to  image  or  detect  tunnels  in  the  presence  of  clutter 
remains.  The  general  problem  of  imaging  in  a random  medium  has  been  studied 
[55],  but  the  specific  application  of  tunnel  detection  might  require  special 
techniques  to  enhance  the  tunnel  signal  and  to  reduce  the  clutter  effects. 
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It  seems  likely  that  a combination  of  further  analysis  coupled  with  EM 
measurements  at  a well  characterized  site  will  yield  the  most  rapid  progress 
on  the  overall  tunnel  detection  problem. 
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APPENDIX  A - VOLUME  INTEGRATION  FOR  ELECTRIC  DIPOLE  SOURCE 


The  unit  vectors  in  (28)  can  be  written  in  terms  of  the  Cartesian  unit 
vectors : 


A A A A 

# = x cos0cos#  + y sin0cos#  - z sin# , 

(Al) 

A A A 

0 = -x  sin0  + y cos 0. 


We  can  write  the  primed  unit  vectors  in  (29)  in  a similar  manner: 


A A A A 

#'  = x cos0'cos#'  + y sin0'cos#'  - z sin#', 

(A2) 

A A A A 

r'  = x cos0'sin#'  + y sin0'sin#'  + z cos#'. 


Using  (Al)  and  (A2) , we  can  evaluate  the  dot  products  required  in  (28) : 


A A 

#•#'  = cos  (0-0 ' ) cos#cos#  ' + sin#sin#', 


A A 

#*r'  = cos  (0-0 ' )cos#sin#  ' - sin#sin#', 

(A3) 

A A 

0«#'  = sin(0  ' -0)  cos#  ' , 


A A 

0*r ' = sin(0 ' -0) . 


The  volume  integral  in  (28)  can  be  written 


A r\  C 2 7T  7T  A 0 0 

/ 1#  • E (r')|  dr'  = f f f |#  • E (r')|  r'  sin#'  d#'d0'dr'.  (A4) 

V ^ b 0 0 " U 

0 0 


The  r'  integration  involves  only  powers  of  r'  and  can  be  done  analytically. 
The  0'  and  #'  integrations  involve  only  trigonometric  functions,  and  they 
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can  also  be  done  analytically.  When  we  substitute  (Al)-(A3)  and  (29)  into 
(A4)  and  carry  out  the  r' , <f> ' , and  9'  integrations,  we  obtain  the  final 
result  in  (30) . 


APPENDIX  B - VOLUME  INTEGRATION  FOR  MAGNETIC  DIPOLE  SOURCE 

The  dot  products  involving  the  unit  vectors  in  (41)  can  be  written 


#•  (rx<£'  ) = <y  • (#xr)  = -<f>' 

(Bl) 

AAA  A AA  AA 

<f>»(rx4>')  = r)  = 4>'»9. 


The  dot  products  in  (Bl)  are  evaluated  in  the  same  manner  as  in  Appendix  A, 
and  the  results  are 


-<f> ' »<f>  = -cos  (<£-<£')  , 


(B2) 


4>'*9  = sin(<f> -<f> ' ) cos#  . 


The  volume  integral  in  (41)  takes  the  same  general  form  as  that  in  (A4) , and 
again  the  r ' , <f>'  , and  6'  integrations  can  be  done  analytically.  When  we 
substitute  (Bl) , (B2) , and  (42)  into  (41)  and  carry  out  the  integrations,  we 
obtain  the  final  result  in  (41) . 


APPENDIX  C - ANGULAR  INTEGRATION  FOR  SCATTERING  CROSS  SECTION 


The  first  step  in  the  evaluation  of  (55)  is  to  write  the  solid  angle 
integration  in  spherical  coordinates: 


2tt  tt 

pa  = f f a sin#  d#  d <f> , (Cl) 

S 0 0 
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where  o is  given  by  (20)  and  the  geometry  is  shown  in  figure  3.  The  sin  x 
factor  can  be  written  in  terms  of  and  9 : 

sin^x  = 1 - cos^  sin^0.  (C2) 

For  small  values  of  k^a,  we  can  approximate  (Cl)  in  the  following  form: 


P° . 


. 4 3 2 _ 

k a <€-,>  2 7T  7T 

J f (1  - cos  <f>  sin  6)  sin0  d 9 d<f> . 


2n 


0 0 


(C3) 


The  9 and  <j>  integrations  in  (C3)  can  be  done  in  closed  form  to  yield  the 
desired  final  result  for  small  k^a: 


4 3 2 

pa  * (4/3)  ktTa  <e7>. 
s 0 1 


(04) 


2 

When  k^a  is  large,  we  make  small-angle  approximations  for  sin  9 and 

2 

sin  (9/2),  and  (Cl)  can  be  approximated 


P°. 


i 4 3 2 0 

kna  <e->  2n  A 

0 j; 


9 d 9 d6 


2tt 


2 2 2 2 

0 0 (1  + k^a  9 ) 


(C5) 


where  A is  chosen  large  enough  to  include  the  significant  portion  of  the 
integrand,  but  small  enough  for  the  small  9 approximations  to  be  valid.  The 
<f>  integration  in  (C5)  contributes  a 2n  factor.  The  9 integration  is 
evaluated  by  substitution: 


A 

; 


9 d 9 


A ,,  .2  2.2.2 
0 (1  + k^a  9 ) 


A , 

p dx 

J 2 2 2 

0 2(1  + k^a  x) 


2 2 
2k0a 


(1  - 


1 + 


v2  2a2 
k^a  A 


) 


(06) 
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We  can  neglect  the  second  term  in  the  parentheses,  and  the  desired  final 


result  for  large  k^a  is 


P° . 


2 2 
k0a  <ei>/2- 
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Clutter 


O 0 


Figure  1.  Typical  geometry  for  tunnel  detection  in  the  presence  of  clutter. 
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Random 


Receiver 


Figure  2.  Geometry  for  scattering  by  an  elemental  volume  dV  in  a random 
medium . 
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A 


Figure  3.  Geometry  for  the  scattering  angle  6 and  the  polarization  angle  \ • 

A A 

The  unit  vectors  i and  o indicate  the  propagation  directions  of 

A 

the  incident  and  scattered  waves,  and  the  unit  vector 
indicates  the  polarization  of  the  incident  wave. 
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Figure  4.  An  electric  dipole  source  at  the  center  of  a random  medium  with 
spherical  boundaries. 
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Figure  5 . Radiation  pattern  of  the  incoherent  scattered  field  for  an 
electric  dipole  source. 
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Figure  6.  Radiation  pattern  of  the  incoherent  field  of  an  electric  dipole 
for  a smaller  value  of  k^b . 
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Figure  7.  Radiation  pattern  of  the  incoherent  field  of  a magnetic 
source . 


dipole 
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k0b 

Figure  8.  Ratio  of  radiated  incoherent  to  coherent  power  for  electric  and 
magnetic  dipole  sources. 
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Figure  9 


Random 

Medium 


Geometry  for  dipole- to-dipole  transmission  in  a random  medium. 
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